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ABSTRACT 

This  paper  presents  a  search  for  radio  transients  at  a  frequency  of  73.8  MHz  (4  nr  wavelength) 
using  the  all-sky  imaging  capabilities  of  the  Long  Wavelength  Demonstrator  Array  (LWDA).  The 
LWDA  was  a  16-dipole  phased  array  telescope,  located  on  the  site  of  the  Very  Large  Array  in 
New  Mexico.  The  field  of  view  of  the  individual  dipoles  was  essentially  the  entire  sky,  and  the 
number  of  dipoles  was  sufficiently  small  that  a  simple  software  correlator  could  be  used  to  make 
all-sky  images.  From  2006  October  to  2007  February,  we  conducted  an  all-sky  transient  search 
program,  acquiring  a  total  of  106  hr  of  data;  the  time  sampling  varied,  being  5  minutes  at  the 
start  of  the  program  and  improving  to  2  minutes  by  the  end  of  the  program.  We  were  able  to 
detect  solar  flares,  and  in  a  special-purpose  mode,  radio  reflections  from  ionized  meteor  trails 
during  the  2006  Leonid  meteor  shower.  We  detected  no  transients  originating  outside  of  the  solar 
system  above  a  flux  density  limit  of  500  Jy,  equivalent  to  a  limit  of  no  more  than  about  10-2 
events  yr-1  deg-2,  having  a  pulse  energy  density  >  1.5  x  10-20  J  nr  2  Hz-1  at  73.8  MHz  for 
pulse  widths  of  about  300  s.  This  event  rate  is  comparable  to  that  determined  from  previous 
all-sky  transient  searches,  but  at  a  lower  frequency  than  most  previous  all-sky  searches.  We 
believe  that  the  LWDA  illustrates  how  an  all-sky  imaging  mode  could  be  a  useful  operational 
model  for  low-frequency  instruments  such  as  the  Low  Frequency  Array,  the  Long  Wavelength 
Array  station,  the  low-frequency  component  of  the  Square  Kilometre  Array,  and  potentially  the 
Lunar  Radio  Array. 
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1.  Introduction 

Transient  emissions—  in  the  form  of  bursts, 
flares,  and  pulses  from  compact  sources — are  the 
signposts  for  explosive  or  dynamic  events.  As 
such,  transient  sources  offer  insight  into  a  variety 
of  fundamental  aspects  of  physics  and  astronomy, 
ranging  from  studying  the  mechanisms  of  particle 
acceleration  on  the  Sun  and  nearby  stars  to  track¬ 
ing  stellar  evolution  and  death  across  the  Universe 
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to  probing  the  intervening  medium(a). 

At  radio  wavelengths,  there  are  well-known 
classes  of  transients,  such  as  the  Sun  and  ra¬ 
dio  pulsars,  as  well  as  a  long  history  of  observ¬ 
ing  transients  from  triggers  at  other  wavelengths, 
such  as  7-ray  burst  (GRES)  afterglows  (e.g., 
Poolev  et  al.  1997:  Tavlor  et  al.  1998:  iFrail  et  ahl 
2000I).  or  the  monitoring  of  known  transient 
sources  such  as  X-ray  binaries  and  microquasars 
'e.g.,  Waltman  fe;  Ghig<3_  19981:  Klein- Wolt  et  al 


2002:  Marti  et  all 20021  Ruuen  et  al.l  20o5b  Fur¬ 
ther,  a  series  of  observations  and  discoveries 
over  the  past  decade  have  emphasized  that 
the  radio  sky  may  be  quite  dynamic.  Known 
sources  have  been  discovered  to  behave  in  new 
ways  and  what  may  be  entirely  new  classes 
of  sources  have  been  discovered — pulsed  ra¬ 
dio  emission  has  been  observed  from  brown 
dwarfs  I  Hallinan  et  al.  20071)  and  formerly  radio¬ 
quiet  magnetars  (Camilo  et  al.l  [2006h;  single  or 
highly  intermittent  pulses  have  been  detected 
from  neutron  stars,  also  known  as  rotating  ra¬ 


dio  transients  (RRATs.  iMcLaughlin  et  al 


Keane  et  al.  201ft  Burke-Soolaor  fc  Bailesl  2010h: 


iiol 


intense  giant  pulses  have  been  detected  from  the 
Crab  pulsar  (Hankins  et  al.  20031):  and  several 
as-yet  unidentified  radio  transients  have  been 
found  (Hvman  et  al.  200.4  Bower  et  al.  2007: 
Lorimer  et  alJ  2007). 

There  are  also  numerous  classes  of  objects  that, 
by  extension  of  known  physics,  have  been  hypoth¬ 
esized  to  be  radio  emitting,  and  which  might  ap¬ 
pear  as  radio  transients.  These  include  extrasolar 
planets  (Farrell  et  al.  1999h.  magnetar  flares  or  gi¬ 
ant  pulses  from  pulsars  in  other  galaxies  (e.g., 
McLaughlin  &:  Corded  20031) .  prompt  emission 


from  GRBs  ( Usov  fc  Kat d  2000:  Sagiv  &;  Waxman 

2002).  evaporating  black  holes  dllees  1977),  and 
extraterrestrial  transmitters  ( Cordes  et  al  •  1 997t). 


Meter- wavelength  observations  have  a  long  his¬ 
tory  of  being  used  for  studying  the  time  domain. 
For  example,  solar  radio  emissions  have  long  been 
known  to  be  dynamic,  and  pulsars  were  first  dis¬ 


covered  in  81  MHz  observations  (jHewish  et  al 
1968h.  Advantages  of  low-frequency  instruments 


for  surveying  the  transient  sky  include  naturally 
large  instantaneous  fields  of  view  and  sensitiv¬ 
ity  to  steep  spectral  index  sources  such  as  might 
result  from  coherent  emission  processes,  while  dis¬ 
advantages  include  radio- wave  propagation  effects 


such  as  dispersion  and  multi-path  propagation, 
absorption  either  within  the  source  or  along  the 
line  of  sight,  ionospheric  disturbances,  and  ra¬ 
dio  frequency  interference  (RFI).  A  number  of 
new  telescopes  with  meter-wavelength  capabili¬ 
ties  are  under  construction,  including  the  Murchi¬ 
son  Wide- field  Array  (MWA,  Lonsdale  et  al 


2009),  the  Low  Frequency  Array  (LOFAR),  the 


Precision  Array  to  Probe  the  E 

poch  of  Reion- 

ization  (PAPER.  Parsons  et  al. 

2010t) .  and  the 

Long  Wavelength  Array  (LWA, 

lEllingson  et  al. 

20091) — many  of  the  science  cases  for  these  tele¬ 


scopes  explicitly  include  searches  for  radio  wave¬ 
length  transients.  It  is  anticipated  that  the  suc¬ 
cess  of  these  telescopes  will  motivate  the  low- 
frequency  component  of  the  Square  Kilometre  Ar¬ 
ray  (SKA,  lDewdney_gt_aL _  200E )  and  the  Lunar 


Radio  Array  (LRA,  lLazio  et  al.  2009t). 


This  paper  presents  an  all-sky  monitoring  cam¬ 
paign  for  radio  wavelength  transients  conducted 
on  the  Long  Wavelength  Demonstrator  Array 
(LWDA).  The  objective  is  two-fold,  both  to  detect 
or  constrain  radio  transients  from  the  observations 
and  to  illustrate  how  future  instruments  could  be 
used  for  all-sky  searching  and  monitoring  obser¬ 
vations.  The  plan  of  this  paper  is  as  follows.  In 
Sj2]  we  describe  the  LWDA  itself  and  its  data  ac¬ 
quisition  path,  in  (J3]we  describe  the  LWDA  data 
acquisition  and  imaging  pipeline  specific  to  the 
formation  of  essentially  all-sky  images  and  radio 
transient  searching,  in  (J3]  we  present  our  results 
and  analysis,  in  fj5]we  discuss  possibilities  for  fu¬ 
ture  instruments,  and  in  Ij6]  we  summarize  our 
conclusions. 


2.  The  Long  Wavelength  Demonstrator 
Array 

The  Long  Wavelength  Demonstrator  Array 
(LWDA)  was  a  testing  platform  developed  by  the 
Naval  Research  Laboratory  prior  to  construction 
of  the  LWA.  The  LWDA  site  work  included  field¬ 
ing  prototype  hardware  for  the  LWA,  development 
of  site  preparation  techniques,  continual  monitor¬ 
ing  of  the  RFI  environment,  software  develop¬ 
ment,  and  initial  science.  Located  on  the  Plains 
of  San  Agustin  in  New  Mexico  near  the  center 
of  the  National  Radio  Astronomy  Observatory’s 
Very  Large  Array  (VLA),  the  LWDA  consisted 
of  16  dual-polarization  dipole  stands  operating 
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Fig.  1.  -  Photo  of  the  Long  Wavelength  Demonstrator  Array  in  New  Mexico  showing  the  individual  dual¬ 
polarization  stands.  Each  stand  consists  of  four  droopy  blades  which  provide  dual-polarization  capabilities. 


as  a  phased  dipole  array  with  a  frequency  range 
of  60-80  MHz  ( York  et  all  20071  Figure  [TJ  ).  The 
element  locations  represented  a  compact  subset 
of  a  larger  256  stand  pseudo-random  distribution 
that  was  developed  for  the  LWA  (IKogan  fc  Cohen 
201)!)..  The  maximum  baseline  for  the  LWDA 
was  approximately  20  m  (Figure  [2]),  providing 
an  angular  resolution  at  the  zenith  of  about  12° 
at  74  MHz. 

Each  of  the  16  LWDA  dipole  stands  consisted 
of  2  “droopy,”  fat  dipoles,  one  per  polarization, 
with  the  dipoles  mounted  at  a  45°  droop  angle  to 
a  central  post  (Figure  [3]).  Each  dipole  consisted 
of  two  “blades,”  each  28  cm  wide  by  1.05  m  long. 
This  design  was  chosen  to  increase  the  symme- 
try  in  the  E-  and  H-planes  (Kerkhoff  fo  Ellingsonl 
20051).  broaden  the  antenna  beam  pattern  (i.e.,  in¬ 


crease  antenna  sky  coverage),  and  obtain  a  larger 
operational  bandwidth. 

The  radio  frequency  (RF)  signals  from  each 


dipole  were  amplified  by  an  active  balun  (24  dB 
of  fixed  gain)  before  being  passed  to  the  LWDA 
electronics,  located  within  a  central,  shielded  shel¬ 
ter.  Inside  the  shelter,  the  incoming  RF  data 
stream  entered  the  digital  receiver  signal  process¬ 
ing  chain.  The  signals  were  digitized  by  a  dual¬ 
channel  10-bit  analog-to-digital  converter  (ADC) 
at  100  megasamples  per  second. 

The  LWDA  electronics  had  the  capability  to 
provide  digital  delay  beamforming  of  two  fully  in¬ 
dependent  beams  of  bandwidth  1.6  MHz  each  as 
well  as  the  capability  for  all-sky  monitoring.  This 
capability  was  implemented  through  two  identi¬ 
cal  signal  processing  chains  that  had  independent 
control  parameters  to  allow  for  two  frequencies, 
two  spatial  beams,  or  two  polarizations.  The  first 
stage  of  digital  processing  was  a  first-in-first-out 
(FIFO)  pipe  that  provided  an  integer  sample  de¬ 
lay  from  0  ns  to  150  ns  in  steps  of  10  ns.  This 
delay  provided  partial  compensation  of  geometric 
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Fig.  2. —  Location  of  the  antennas  within  the 
LWDA. 


path  and  system  time  delays  at  this  stage.  Fol¬ 
lowing  the  FIFO  was  a  complex  mixer  where  the 
in-phase  (I)  and  quadrature  (Q)  components  for 
each  signal  chain  were  generated  using  a  CORDIC 
rotator.  The  signals  then  passed  through  a  low- 
pass  filtering  stage  using  a  Cascaded  Integrator 
Comb  (CIC)  filter  and  were  decimated  by  a  fac¬ 
tor  of  14.  Next,  the  signals  passed  through  an¬ 
other  FIFO,  which  provided  coarse  sample  delays 
in  excess  of  2  /rs  in  140  ns  increments.  The  final 

1.6  MHz  bandpass  for  each  data  stream  was  set 
in  a  finite  impulse  response  (FIR)  filter  that  also 
reduced  the  data  rate  by  another  factor  of  3.  Fol¬ 
lowing  this  stage,  the  data  were  interleaved  and 
passed  to  an  adder  board  at  the  final  data  rate  of 
approximately  4.8  MSPS.  The  adder  board  was 
used  to  either  sum  (beamform)  or  interleave  (all¬ 
sky  image)  the  incoming  coherent  signals. 

The  LWDA  frequency  range  was  chosen,  in 
large  part,  because  early  plans  for  the  array  in¬ 
cluded  the  goal  of  operating  it  as  one  of  the  ele¬ 
ments  of  the  VLA  74  MHz  system  (Kassim  et  al.l 
.2007).  which  in  turn  was  centered  on  the  (primary) 
frequency  allocation  for  radio  astronomy  at  73.0- 

74.6  MHz  in  the  U.S.  The  operational  frequency 
range  was  much  broader  and  includes  frequency 
allocations  for  a  number  of  different  services,  pri¬ 
marily  TV  broadcasting,  at  the  time  of  the  LWDA 
operation  0  Tests  with  the  LWDA  and  the  pro- 

1  Since  the  time  of  these  experiments,  the  U.S.  has  converted 


Fig.  3. —  Photo  of  several  LWDA  elements  show¬ 
ing  the  droopy  blades  and  an  NRAO  VLA  antenna 
in  the  background. 


totype  LWA  equipment  have  shown  the  presence 
of  TV  signals,  but  few  other  strong  emissions  in 


this  band  (Dufhn  &  Rav  2007: 

Hicks  et  al. 

2007: 

Paravastu  et  al.  20071  Ray  et  al. 

20071:  Crai  si 

2008; 

Jaeger  &  Mutel  not  dated).  Further,  experience 

with  the  VLA  74  MHz  system  indicated  that,  for 
that  system,  the  most  significant  source  of  inter¬ 
ference  was  the  VLA  itself.  Indeed,  in  comparison 
with  the  VLA,  for  which  the  74  MHz  feeds  were 
mounted  near  prime  focus  (several  meters  above 
the  ground),  the  LWDA  should  be  less  susceptible 
to  RFI.  Our  experience  with  the  LWDA  was  con¬ 
sistent  with  the  experience  from  the  VLA,  in  that 
no  significant  amount  of  observing  time  was  lost 
due  to  RFI.  (See  also  H4.ll") 


The  LWDA  system  was  sky-noise  dominated  by 
at  least  6  dB  over  its  full  operating  range.  The 
LWDA  construction  was  completed  and  first-light 
images  were  obtained  on  2006  October  23. 

Our  focus  here  is  on  the  all-sky  observations 
obtained  by  the  LWDA,  but  it  also  had  a  two- 
element  interferometer  mode.  In  this  mode,  an 
outlier  dipole  could  be  utilized  to  monitor  the  flux 
densities  of  isolated  point  sources,  contributing 
to  a  study  of  the  secular  flux  density  decrease  of 
Cas  A  (Helmboldt  fe  Kassimll2009l). 


from  analog  to  digital  TV  broadcasts,  and  much  of  the 
LWDA  operating  frequency  range  is  no  longer  used  for  TV 
broadcasts  at  the  time  of  writing  (|Cranell2009h. 


4 


3.  All-Sky  Observations 


The  observations  presented  here  consist  of  op¬ 
erating  the  LWDA  in  its  “all-sky”  imaging  mode 
observing  at  zenith.  In  this  operational  mode, 
the  all-sky  data  acquisition  system  cycles  serially 
through  all  120  pairs  of  dipoles  in  the  array  and 
correlates  the  data  from  each  pair.  Each  base¬ 
line  is  measured  for  51  ms  with  a  total  cycle  time 
of  13  s  for  both  polarizations  of  all  baselines  (in¬ 
cluding  auto-correlations).  As  a  result  of  various 
processing  overheads,  the  on-sky  cadence  of  imag¬ 
ing  described  herein  is  slower.  Our  observations 
began  with  a  5  minute  cadence  (i.e. ,  13  s  to  acquire 
the  full  set  of  cross  correlations  to  form  an  image 
with  images  repeated  every  5  minutes),  which  was 
later  improved  to  a  2  minute  cadence  for  most  of 
the  data  reported  here. 


Our  focus  on  an  (all-sky)  imaging  transient 
search  stems  from  the  fact  that  the  number  of 
dipoles  in  the  LWDA  was  relatively  small  so  that 
cross-correlation  of  all  of  them  was  computation¬ 
ally  tractable.  The  alternative  to  an  imaging 
search,  in  which  transients  are  identified  by  com¬ 
paring  images  of  the  sky  at  different  epochs,  is  a 
non-imaging  search,  in  which  transients  are  iden¬ 
tified  in  time  series  of  the  voltages  or  intensities 
from  a  telescope(s).  Traditionally,  the  tradeoffs 
between  the  two  kinds  of  searches  have  been  that 
imaging  searches  obtain  higher  angular  resolution 
at  the  cost  of  time  resolution  with  the  converse 
being  true  for  non-imaging  searches.  The  choice 
between  the  kind  of  transient  search  to  employ 
depends  not  only  on  the  instrument  available,  but 
also  upon  what  is  known  about  the  transient  pop¬ 
ulation^)  of  interest  ({Nemiroff  1200311. 


The  observations  were  centered  at  a  frequency 
of  73.8  MHz  with  80  spectral  sub-bands  over  a  to¬ 
tal  bandwidth  of  1.6  MHz  (20  kHz  per  sub-band). 
During  the  course  of  the  array  installation,  the 
amplitude  and  phase  of  the  signal  path  from  each 
dipole  were  determined.  These  values  were  applied 
in  real  time  as  the  cross-correlations  (visibilities) 
were  formed.  The  normalized  bandpass  shapes 
were  applied  to  correct  for  individual  telescope 
variations.  The  visibilities  were  then  recorded  on 
disk  for  later  processing. 

Post-processing  consisted  of  the  following  steps: 


1.  The  raw  visibilities  were  converted  into  the 


FITS-IDI  format  dClarke  et  al.ll2007h. 


2.  The  visibility  files  were  ingested  by  the 
Astronomical  Image  Processing  System 
(AIPS)H 

3.  Standard  coplanar  imaging  procedures  were 
used  to  convert  the  visibilities  to  images. 

4.  Strong  sources  in  the  images  were  blanked 
(see  below),  and,  because  the  images  for¬ 
mally  extend  beyond  the  horizon,  they  were 
clipped  at  the  horizon. 


Figure  0]  presents  a  typical  all-sky  image.  Sources 
easily  detected  in  the  LWDA  all-sky  maps  include 
Cas  A,  Cyg  A,  the  inner  Galactic  plane,  North 
Polar  Spur,  Loop  III,  and  the  Sun. 


We  used  the  apparent  brightness  of  both  Cas  A 
and  Cyg  A  to  obtain  a  crude  estimate  of  the  power 
pattern  of  the  LWDA.  Simulations  of_the  power 
pattern  of  the  individual  dipoles  ([York  et  al.]|2007h 
suggest  that  the  gain  decreases  toward  the  horizon 
(as  would  be  expected  from  a  dipole)  as  well  as  to¬ 
ward  the  zenith.  The  latter  is  caused  by  the  design 
goal  of  operating  the  dipoles  over  a  relatively  large 
frequency  range.  The  brightness  measurements  of 
Cas  A  and  Cyg  A  do  show  these  large-scale  gain 
variations  as  well  as  smaller  variations,  which  we 
estimate  to  be  approximately  a  factor  of  2  (3  dB) 
in  amplitude. 


The  all-sky  transient  data  were  passed  into  an 
analysis  pipeline  written  within  the  AIPS  software 
system.  The  pipeline  first  clipped  high  phase- 
center  amplitudes  that  are  known  to  be  generated 
by  sporodic  monitor  and  control  (M&C)  software 
issues.  The  data  were  then  Fourier  transformed 
into  dirty  images  of  the  entire  sky  visible  at  the 
LWDA  site.  Even  though  the  LWDA  can  ob¬ 
serve  essentially  the  entire  sky,  “wide-field”  imag¬ 
ing  techniques  ( Cornwell  fc  Perlevl  1992)  are  not 
required  given  the  size  of  the  LWDA,  its  opera¬ 
tional  frequency,  and  our  use  of  snapshot  imaging. 


Images  were  blanked  around  the  positions  of 
the  strong  sources  that  could  be  detected  easily 
(i.e.,  Cyg  A,  Cas  A,  the  Sun,  the  inner  Galac¬ 
tic  plane,  Loop  III,  and  the  North  Polar  Spur). 
This  strategy  of  blanking  the  images  allowed  us 
to  use  a  relatively  simple  threshold  test  to  search 
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Fig.  4. —  Comparison  of  all-sky  LWDA  images  (green  contours)  with  the  408  MHz  all-sky  image  (color)  from 


Haslam  et  al.l  (119821).  Insets  in  both  panels  (demarcated  by  the  white  circle)  show  the  all-sky  LWDA  images 


(in  color);  in  these  insets,  north  to  up  and  east  is  to  the  left.  The  dotted  white  line  indicates  the  horizon  for 
the  LWDA,  with  the  portion  above  the  dotted  white  line  being  visible.  The  LWDA  images  are  at  73.8  MHz 
and  were  acquired  on  2006  October  28.  Top:  a  snapshot  acquired  when  the  Galactic  plane  was  passing 
nearly  overhead.  Individual,  well  known  sources  and  the  inner  Galactic  plane  dominate  the  image.  Bottom: 
a  snapshot  acquired  when  the  inner  Galactic  plane  was  below  the  horizon.  While  the  LWDA  image  would 
appear  to  be  only  noise,  there  are  features  that  can  be  identified  with  the  North  Polar  Spur  and  Loop  III 
radio  structures. 
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for  the  presence  of  transients,  at  the  cost  of  reduc¬ 
ing  somewhat  our  sky  coverage.  Near  the  end  of 
the  LWDA’s  operation,  we  did  attempt  to  change 
the  array  so  that  the  time  sampling  was  in  sidereal 
time.  Had  this  been  successful,  it  would  then  have 
been  possible  to  subtract  images  acquired  on  dif¬ 
ferent  days  and  construct  “difference”  images  for 
which  this  blanking  scheme  probably  would  not 
have  been  necessary. 

As  a  measure  of  our  ability  to  detect  a  tran¬ 
sient,  particularly  in  light  of  our  non-simultaneous 
sampling  of  the  baselines,  we  conducted  a  num¬ 
ber  of  simulations  in  which  transients  of  differ¬ 
ent  strengths  and  positions  within  the  image  were 
inserted  into  the  data.  Table  |T|  summarizes  the 
results:  Column  1  indicates  what  fraction  of  the 
baselines  are  illuminated  by  the  simulated  tran¬ 
sient,  Column  2  shows  the  implied  duration  of 
the  transient,  and  Column  3  shows  the  approx¬ 
imate  signal-to-noise  ratio  (S/N)  that  the  tran¬ 
sient  would  have  to  have  in  order  to  be  detected. 
One  important  aspect  of  these  simulations  is  that 
the  simulated  transients  were  inserted  into  the  vis¬ 
ibility  data.  As  such,  these  simulations  do  not 
take  into  account  the  antenna  power  pattern.  The 
trend  in  Table  [1]  of  increasing  S/N  required  in  or¬ 
der  to  detect  the  transient  is  perhaps  not  surpris¬ 
ing,  but  we  find  it  reassuring  that  we  could  detect 
transients,  even  those  whose  durations  are  only  a 
minor  fraction  of  the  time  it  took  the  LWDA  to 
acquire  a  full  sample  of  visibility  data. 

Using  our  simulations  and  visual  inspection  of 
the  images  as  a  guide,  each  all-sky  snapshot  image 
was  searched  for  events  above  the  5er  level.  Any 
image  containing  a  potential  transient  event  was 
saved  for  further  analysis. 

4.  Results 

We  begin  this  section  by  both  motivating  and 
illustrating  the  capability  of  the  LWDA  to  detect 
radio  transients,  and  we  then  turn  to  the  detection 
of  astronomical  radio  transients. 


4.1.  Detected  Radio  Transients:  The  Sun 
and  Meteor  Trails 


The  Sun  is  a  well-known  variable  radio  source 
at  tens  of  MHz  (e.g.,  IWild  fe  McCreadv  1950lf. 
While  the  Sun  is  currently  in  a  state  of  excep¬ 
tionally  low  radio  emission,  during  the  course  of 


the  LWDA  observations  reported  here,  a  number 
of  solar  radio  bursts  were  detected,  including  a  se¬ 
ries  of  intense  radio  bursts  on  2006  December  14 
(Figure  [5]),  probably  associated  with  the  X-class 
flares  that  occurred  then. 

The  focus  of  our  monitoring  program  was  pri¬ 
marily  on  possible  extrasolar  system  transients, 
and,  as  such,  we  made  no  effort  to  track  the  solar 
flux  density  (32]).  However,  a  similar  monitoring 
program  could  be  used  to  track  the  flux  density 
of  the  Sun  at  these  frequencies.  Such  a  program 
would  be  more  valuable  if  the  array  monitoring 
the  Sun  had  both  a  larger  bandwidth  and  a  higher 
time  resolution.  With  these  improved  characteris¬ 
tics,  not  only  would  such  an  array  be  capable  of 
detecting  solar  radio  bursts,  it  might  be  able  to 
image  their  evolution  in  frequency. 

Conversely,  as  Figure  [5]  illustrates,  the  Sun  can 
become  a  source  of  “interference.”  During  these 
intense  solar  radio  bursts,  the  dynamic  range  of 
the  LWDA  and  imaging  pipeline  were  sufficiently 
limited  that  only  the  Sun  was  visible.  Whether 
other  observations  can  be  conducted  during  so¬ 
lar  radio  bursts  will  depend  upon  the  imaging  dy¬ 
namic  range  of  future  systems,  however,  exceeding 
an  imaging  dynamic  range  of  order  100  seems  dif¬ 
ficult  for  rapid,  all-sky  imaging  pipelines. 

Another  well-known  class  of  radio  transients 
is  reflections  from  ionized  meteor  trails.  As  a 
meteor  enters  Earth’s  atmosphere,  its  velocity  is 
high  enough  that  it  can  produce  an  ionized  trail, 
which  in  turn  can  reflect  a  radio  transmission  (e.g., 
iMillman  et  al.l[l948). 

In  2006  November,  during  the  Leonid  meteor 
shower,  data  were  acquired  with  the  LWDA  tuned 
to  a  frequency  of  61  MHz.  At  the  time  of  these  ob¬ 
servations,  this  frequency  was  within  the  band  al¬ 
located  for  television  broadcasting,  and  there  were 
a  number  of  TV  stations  in  New  Mexico  and  the 
southwest  U.S.  that  used  this  allocation.  It  is  a 
simple  matter  to  show  that  even  a  TV  transmit¬ 
ter  of  relatively  modest  power  located  more  than 
1000  km  distant  would  still  produce  a  strong  re¬ 
flection.  For  the  purposes  of  this  observation,  the 
process  for  data  acquisition  was  modified  to  re¬ 
duce  the  number  of  bits  and  allow  simultaneous 
sampling  of  all  16  dipoles.  Figure  [6]  illustrates 
two  examples  of  a  meteor  reflection  imaged  by  the 
LWDA. 
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Table  1 

Transient  Detection  Simulations 


Baseline  Fraction 
Illuminated 

Time 

(s) 

S/N 

100% 

6 

5 

90% 

5.4 

7 

75% 

4.5 

10 

50% 

3 

13 

30% 

1.8 

30 

Note. — The  results  of  injecting 
simulated  transients  into  the  LWDA 
data.  Column  1  indicates  what  frac¬ 
tion  of  the  baselines  are  illuminated 
by  the  simulated  transient.  Because 
the  LWDA  acquired  all-sky  imag¬ 
ing  data  by  sampling  baselines  se¬ 
quentially,  Column  2  shows  the  im¬ 
plied  duration  of  the  transient.  Col¬ 
umn  3  shows  the  approximate  signal- 
to-noise  ratio  (S/N)  that  the  transient 
would  have  to  have  in  order  to  be  de¬ 
tected. 


Fig.  5.-  Illustration  of  the  LWDA’s  detection  of  a  solar  radio  burst.  Left:  the  sky  before  the  solar  radio 
burst.  Right:  the  sky  during  the  solar  radio  burst  of  2006  December  14.  This  example  is  one  of  the  several 
intense  radio  bursts  that  occurred  in  2006  December. 
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Fig.  6.  LWDA  images  of  reflected  TV  signals  from  ionized  meteor  trails  during  the  2006  Leonid  meteor 
shower.  Each  image  is  constructed  from  a  0.1  s  data  acquisition  at  a  frequency  of  61  MHz.  Top  left  and  top 
right:  images  of  reflected  signals  from  Leonid  meteors.  Bottom:  a  reference  image  showing  the  nominal  sky 
observed  by  the  LWDA.  The  zenith  is  at  the  center  of  the  image.  The  two  discrete  emitters  near  the  center 
of  the  field  are  celestial  radio  sources  (Cas  A  and  Cyg  A).  The  band  stretching  from  the  upper  right  to  lower 
left,  terminating  in  the  strong  emission  near  the  edge  of  the  field  is  from  the  Milky  Way. 
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4.2.  Astronomical  Radio  Transients 


10  kpc 


1  Mpc 


The  observing  campaign  lasted  from  2006  Oc¬ 
tober  27  to  2007  February  17,  during  which  there 
were  59  days  on  which  some  useful  data  were  ac¬ 
quired.  On  2006  December  13,  we  were  able  to 
improve  the  cadence  at  which  images  were  formed 
from  5  minutes  to  2  minutes.  The  resulting  to¬ 
tal  observing  time  for  the  transient  campaign  was 
106  hr.  During  these  observations,  no  transients 
were  detected  above  the  5cr  level  in  any  image. 
Based  on  comparative  analysis  of  the  brightnesses 
of  Cas  A  and  the  Sun,  we  estimate  the  (la)  noise 
level  in  the  LWDA  all-sky  snapshots  to  be  roughly 
500  Jy.  We  also  assume  that  the  effective  sky 
coverage  of  the  LWDA  is  roughly  10,000  deg2 
(~  7 r  sr).  Although,  we  show  essentially  images 
of  the  entire  hemisphere,  the  gain  of  the  LWDA 
was  unlikely  to  be  constant  across  the  sky,  an  as¬ 
sumption  partially  supported  by  modeling.  Conse¬ 
quently,  we  assume  that  the  effective  sky  coverage 
is  somewhat  less  than  a  full  hemisphere. 


A  total  of  29,437  data  sets  were  processed 
through  the  imaging  pipeline.  Of  these,  1764  pro¬ 
duced  potential  transient  candidates.  Visual  in¬ 
spection  of  these  showed  that  they  could  be  ex¬ 
plained  as  non-astronomical  transients:  slight  er¬ 
rors  in  our  calculations  of  the  extents  of  strong 
sources  resulted  in  one  pixel  of  a  strong  source 
not  being  blanked  appropriately. 


Combining  these  results,  we  place  a  limit  of 
<  10-2  events  yr-1  deg-2,  having  a  pulse  energy 
density  >  1.5  x  10-20  J  m-  Hz-1  at  73.8  MHz 
for  pulse  widths  of  about  300  s.  We  now  con¬ 
sider  multiple  characterizations  of  this  result  for 
easy  comparison  with  existing  results  in  the  liter- 
ature.  Our  focus  is  on  comparison  with  MOTOR 
( Amv  et  al.  1989tl  and  STARE  ( Katz  et  all  20031 . 


as  both  represent  nearly  all-sky  surveys  at  frequen¬ 
cies  below  1  GHz. 

Following  iKatz  et  al  (2003[).  the  LWDA  flux 
density  limit  can  be  recast  in  terms  of  the  bright¬ 
ness  temperature  that  an  object  would  have  to 
have  in  order  to  be  detectable.  Using  parameters 
appropriate  for  the  LWDA,  and  assuming  a  quite 
conservative  5000  Jy  flux  density  limit  (10cr),  we 
find 

Tb>  40  K  ,  (1) 


for  an  object  of  linear  size  L  located  at  a  dis- 


1  pc 


Log[Distance  (cm)] 


Fig.  7.  Brightness  temperature  limits  for  ob¬ 
jects  of  varying  sizes  (diagonal  lines)  as  a  function 
of  the  distance  to  the  object.  Lines  for  fiducial 
object  dimensions  are  shown  as  are  fiducial  dis¬ 
tances.  The  required  brightness  temperature  for 
an  object  to  be  detectable  is  above  and  to  the  left 
of  a  line  of  constant  object  dimension.  The  range 
of  the  ordinate  is  approximately  that  allowed  for 
the  LWDA  observations  (viz.  Figure  0. 


tance  D.  Figure  [T]  shows  the  resulting  bright¬ 
ness  temperature  limits  for  transients  detectable 
by  the  LWDA.  Sufficiently  strong  bursts  from  a 
neutron  star  or  a  “super-flare”  from  a  nearby  star 
are  among  the  possible  sources  that  might  have 
been  detectable  by  the  LWDA  (or  an  instrument 
like  the  LWDA).  We  consider  these  possibilities 
further  in  fj5] 

Table  m  compares  our  rate  with  those  of  other 
“all-sky”  surveys  at  other  wavelengths  from  the 
literature.  In  constructing  Table  [2]  we  have 
restricted  our  focus  to  other,  largely  similar 
search  programs,  namely  those  at  frequencies  be¬ 
low  1  GHz  using  dipole-based  arrays.  We  have 
also  restricted  Table  0  to  programs  that  conducted 
blind  surveys,  as  opposed  to  those  that  targeted 
known  sources  (e.g.,  pulsars,  GRBs). 

In  general,  our  LWDA  search  program  produces 
an  upper  limit  to  the  event  rate  that  is  compara¬ 
ble  to  or  competitive  with  other  search  programs, 
typically  at  higher  frequencies.  However,  our  se¬ 
rial  data  acquisition  scheme  results  in  a  somewhat 
poorer  limit  on  the  pulse  energy  density.  If  the 
LWDA  data  acquisition  system  had  had  a  more 
rapid  sampling  time,  e.g.,  on  a  scale  of  order  0.1  s 
(the  time  scale  for  the  acquisition  of  a  signal  from 
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Table  2 

Transient  Event  Rates 


Program 

Rate  Limit 
(events  yr-1  deg-2) 

Pulse 

Energy  Density 
(J  m-2  Hz-1) 

Timescale 

00 

Frequency 

(MHz) 

Scaled 

Pulse  Energy  Density 
(J  m-2  Hz-1) 

Reference 

LWDA 

10-2 

1.5  x  10-20 

300 

73.8 

this  work 

HM74 

0.05 

~  10-23 

-  0.1 

270 

3.4  x  10-23 

1 

K+77 

4  x  10-3 

2.5  x  10-22 

>  1 

370-550 

1.3  x  10-21 

2 

4  x  10-3 

3  x  10-22 

>  1 

38-60 

1.5  x  10-22 

2 

MOTOR 

163 

Iq-28 

0.025 

843 

1.1  x  10-27 

3 

STARE 

10-3 

10-24 

0.125 

611 

8.3  x  10-24 

4 

Note. — The  scaled  pulse  energy  density  is  the  pulse  energy  density  scaled  to  a  frequency  of  73.8  MHz  assuming  a 
nominal  spectral  index  of  —  1.  References:  (1)  iHuguenin  fc  Moorel  |l974]);  (2)  iKardashev  et  all  <119771^ :  (3)  lAmv  et,  alJ 
(1 198911:  141  iKatz  et  all  ll2003h. 


an  individual  dipole),  our  pulse  energy  density 
limit  would  have  been  comparable  to  many  of  the 
other  search  programs. 

An  alternate  means  of  characterizing  radio 
transient  searches  is  in  terms  of  a  “phase  space” 
diagram  (Cordes  et  al.ll2004blh  which  can  be  cast 
in  terms  of  the  brightness  temperature  of  a  radio 
transient  in  the  Rayleigh- Jeans  limit  (Equation!]]) . 
Figure[5]shows  the  region  of  phase  space  probed  by 
the  LWDA  in  comparison  with  the  range  defined 
by  various  known  and  hypothesized  phenomena. 
In  producing  this  plot,  we  have  assumed  a  lower 
distance  limit  to  any  celestial  radio  transient  to 
be  1  pc,  such  as  might  be  expected  from  a  nearby 
low-mass  star,  brown  dwarf,  or  planet  ((J5]) .  In 
practice,  the  local  density  of  sources  may  be  suffi¬ 
ciently  low  that  there  are  no  transients  within  1  pc, 
which  would  have  the  effect  of  increasing  the  min¬ 
imum  value  of  SD 2  that  the  LWDA  could  have 
probed.  Clearly  the  upper  distance  limit  could  be 
the  edge  of  the  observable  Universe,  which  would 
imply  correspondingly  higher  brightness  tempera¬ 
tures. 


As  a  final  approach  to  constraining  the  event 
rate  for  transients,  we  consider  a  probabilistic  one. 
The  typical  determination  of  an  event  rate  (such 
as  those  in  Table  is  derived  from  the  total  ob¬ 
serving  duration  T  and  the  solid  angle  Q  observed. 
Instead,  we  ask  what  is  the  probability  p  of  detect¬ 
ing  a  transient  in  a  given  observation,  given  that 


no  transients  were  detected  in  our  series  of  obser¬ 
vations. 


We  conducted  29,437  trial  observations  dur¬ 
ing  jthecourse  of  this  observational  program.  As 
DeGroot  ( 19861)  discusses,  the  probability  of  ob¬ 


taining  i  detections  in  N  trials,  given  that  the  de¬ 
tection  probability  in  any  individual  trial  is  p ,  can 
be  related  to  the  beta  distribution 


P(n,j) 


T(n  +  2) 

T(j  +  l)r(n  — j  +  1) 


P°(l  -P)n~j, 


(2) 


where  T(a;)  is  the  gamma  function,  j  =  Xi, 

and  Xj  =  1  for  a  detection  or  0  for  a  non-detection. 
That  is,  j  is  simply  the  count  of  the  total  number 
of  detections. 

With  j  =  0,  and  integrating  Equation  ([2])  to 
form  the  incomplete  beta  function,  we  can  deter¬ 
mine  the  maximum  detection  probability  p  for  a 
specified  confidence  level.  For  the  LWDA  or  an 
array  with  LWDA-like  parameters  («  10,000  deg2 
sky  coverage,  73.8  MHz  frequency,  <  300  s  du¬ 
ration),  we  find  that  the  probability  of  detecting 
a  transient  in  a  single  LWDA  observation  cannot 
exceed  p  =  2  x  10-4  lest  the  probability  that  we 
should  have  detected  a  transient  in  the  full  LWDA 
transient  campaign  exceeds  99.7%. 

Given  the  relatively  large  solid  angle  accessed 
by  the  LWDA,  only  a  modest  improvement  in 
the  probability  of  detecting  transients  would  re¬ 
sult  from  increasing  the  sky  coverage  for  future  in- 
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log  vW  (GHz  s) 

Fig.  8. —  Phase  space  defined  by  radio  transients.  The  abscissa  is  an  uncertainty-like  quantity  given  by  the 
product  of  the  observing  frequency  v  and  the  duration  of  the  transient  W .  For  our  LWDA  survey,  we  define 
two  regions.  The  first  is  from  2  to  6  s  and  is  the  interval  required  to  accumulate  the  visibility  data  from  all  of 
the  interferometric  baselines  (10  and  Table  [lj.  The  second  is  from  5  minutes  to  24  hr  and  covers  the  cadence 
at  which  individual  images  were  made  to  the  longest  continuous  observing  duration  during  our  program. 
The  ordinate  is  a  pseudo-luminosity  Sp^D2,  for  a  peak  flux  density  S^k  and  distance  D.  We  take  a  lower 
distance  limit  of  1  pc,  such  as  might  be  expected  from  a  nearby  low-mass  star,  brown  dwarf,  or  extrasolar 
planet.  Various  classes  of  known  or  hypothesized  transients  are  indicated.  The  sloping  solid  lines  indicate 
lines  of  constant  brightness  temperature,  with  a  brightness  temperature  of  1012  K  taken  as  the  dividing  line 
between  coherent  and  incoherent  emission. 
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struments  (or  would  have  resulted  from  increasing 
the  sky  coverage  of  the  LWDA  itself).  Rather,  an 
increased  probability  is  likely  to  result  only  from 
improved  sensitivity,  developing  a  search  method 
that  does  not  require  blanking  the  Galactic  plane, 
or  both. 


5.  Future  Possibilities 


In  the  near  future,  the  first  station  of  the  Long 
Wavelength  Array  (LWA-1),  the  LOFAR,  and  the 
International  LOFAR  Telescope  (ILT)  will  become 
operational,  with  LOFAR  being  a  subset  of  the 
ILT.  All  will  consist  of  dipole  stations  having  an 
operational  frequency  range  overlapping  or  com¬ 
parable  to  that  of  the  LWDA. 


The  individual  dipole  stations  of  LOFAR  and 
the  ILT  will  have  a  factor  of  a  few  more  dipoles 
than  the  LWDA  while  LWA-1  will  256  dipoles 
(16 x  as  many),  but  they  may  not  be  able  to  im¬ 
prove  upon  our  flux  density  limits  substantially 
because  the  diameters  of  the  LWDA,  LWA-1,  and 
LOFAR  stations  are  sufficiently  small  that  they 
are  confusion  limited  on  short  time  scales.  How¬ 
ever,  with  multiple  LOFAR  stations,  the  identifi¬ 
cation  and  excision  of  RFI  may  be  improved  sub¬ 
stantially  (Bhat  et  al.  2005h .  In  the  case  of  the 


LWA-1,  all-sky  imaging  will  be  carried  out  by  pass¬ 
ing  a  narrow  bandwidth  (~  100  kHz)  signal  from 
each  of  the  256  stands  to  a  dedicated  cluster  for 
correlation  and  imaging.  This  observing  mode  op¬ 
erates  in  tandem  with  the  beamforming  modes  so 
that  it  will  be  possible  to  image  the  visible  hemi¬ 
sphere  full  time. 


The  MWA  will  also  conduct  transient  searches, 
at  a  somewhat  higher  frequency  range  than  the 
LWDA,  but  it  will  likely  not  include  an  all-sky 
imaging  capability.  Building  upon  these  next- 
generation  low-frequency  arrays  will  be  the  SKA 
and  the  LRA.  Their  designs  will  be  influenced  by 
the  work  on  LOFAR,  MWA,  and  similar  low  ra¬ 
dio  frequency  interferometers,  but,  in  both  cases, 
we  anticipate  that  transient  searches,  and  poten¬ 
tially  all-sky  imaging,  will  be  a  capability.  As  Fig¬ 
ure  [H]  illustrates,  and  as  we  now  discuss  below,  in 
principle,  an  imaging  dipole  station,  similar  to  the 
LWDA,  could  probe  a  variety  of  known  or  hypoth¬ 
esized  classes  of  transients. 


5.1.  Radio  Pulses  from  Ultra-high  Energy 
Cosmic  Rays  or  Neutrinos 

Intense,  short-duration  pulses  (~  1  MJy  in 
~  10  ns)  at  decimeter  and  meter  wavelengths 
have  been  detected  from  the  impact  of  ultra-high 
energy  cosmic  rays  on  Earth’s  atmosphere  (e.g., 
Jelley  et  al.  19651:  Weekesl[200li:  Falcke  fc  Gorhaml 

20031).  and  in  at  least  one  case  even  imaged 
by  a  radio  array  not  dissimilar  to  the  LWDA 


( Falcke  et  al.  20051).  High-energy  neutrinos  im¬ 


pacting  the  lunar  regolith  should  also  produce  ra¬ 
dio  pulses  via  the  Askar ’van  effect  (jDagkesamanskii  &  Zheleznvkh 


1989;  iFalcke  fc  Gorhai— 
pulses  have  been  detected 


tr  yan  ertec 
3  l2Q03h. 


though  no  such 
from  the  Moon  to 


date 

(Hankins  et  al. 

1991 

Gorham  et  al.  2004 

Beresnyak  et  al.  2005t 

Langston  et  al. 

2009 

Scholten  et  al. 

2009 

Buitink  et  al.  2010::  James  et  al. 

2010 

Jaeger  et  al. 

20101).  In  the  case  of  the  lunar  neutrino  ra¬ 
dio  pulses,  the  expected  pulse  amplitude  is  more 
poorly  constrained  because  it  depends,  in  part, 
upon  properties  (primarily  surface  roughness 
of  the  lunar  surface 
Scholten  et  alJ  l2006t 


Alvarez-Muniz  et  al 


James  fc  Protheroe 


200' 

20091). 


Nonetheless,  if  sufficiently  energetic  neutrinos  ex¬ 
ist,  it  is  plausible  that  radio  pulses  would  be  emit¬ 
ted  at  frequencies  relevant  for  an  LWDA-like  in¬ 
strument. 

The  intrinsic  width  of  these  radio  pulses  is  ex¬ 
tremely  short,  far  shorter  than  the  51  ms  data 
acquisition  time  for  an  individual  baseline  of  the 
LWDA;  as  an  example,  Falcke  et  al.  ( 20051)  mea¬ 


sure  a  width  narrower  than  30  ns.  As  noted  above 
(Jf3j).  during  the  commissioning  of  the  LWDA,  a 
sporadic  M&C  issue  was  recognized  in  that  occa¬ 
sionally  a  baseline  would  produce  an  anomalously 
large  amplitude.  Upon  inspection,  we  found  that 
baselines  with  anomalously  large  amplitudes  also 
had  a  phase  of  0°  (which  is  one  of  the  reasons 
that  we  explain  these  anomalous  amplitudes  as 
being  due  to  an  M&C  issue).  Further,  unlike  the 
LOFAR  Prototype  Station  (LOPES,  [Falcke  et  al 


120051).  we  had  no  co-located  particle  detector  to 
use  as  a  trigger  for  cosmic  ray  detection.  Subse¬ 
quently,  we  employed  a  simple  threshold  test,  in 
which  baselines  with  anomalous  amplitudes  and 
phases  near  zero  were  excised.  A  consequence  of 
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3  No  pulses  from  the  Moon  have  been  detected,  but  the 
Askar ’yan  effect  has  been  demonstrated  in  terrestrial  ac¬ 
celerators  dGorham  et  al.l[200f)l.  [20071). 


this  flagging  methodology  is  that  radio  pulses  from 
high-energy  particles  at  the  zenith  would  not  be 
recognized  by  our  analysis  of  LWDA  data. 

5.2.  Stars  and  Substellar  Objects 

The  Sun  generates  intense  emission  at  LWDA 
frequencies,  notably  Type  III  and  IV  radio  bursts, 
but  even  the  strongest  such  radio  bursts  are  far  too 
faint  to  be  detected  over  interstellar  distances  with 
the  current  most  sensitive  meter-  or  decameter- 
wavelength  telescopes.  Gereelv  ( 19861)  considered 
the  detection  of  solar-type  stars  at  low  radio  fre¬ 
quencies  and  finds  that  millijansky  sensitivities 
would  likely  be  required  to  detect  the  equiva¬ 
lent  of  the  most  intense  solar  radio  bursts  from 
nearby  stars.  There  are  approximately  20  solar- 
type  stars  in  the  solar  neighborhood  (<  10  pc, 
Holmberg  et  al.l  120031.  The  absence  of  any  tran¬ 
sients  that  could  be  identified  with  these  stars  puts 
only  modest  constraints  on  their  levels  of  activ¬ 
ity.  Not  only  is  the  LWDA  sensitivity  likely  to 
be  insufficient,  if  these  stars  have  solar-like  cycles 
over  decadal  time  scales,  the  absence  of  any  flar¬ 
ing  could  merely  reflect  that  a  star  is  currently  in 
a  quiescent  phase  of  the  cycle  (much  like  the  Sun 
is  in  the  current  solar  cycle). 

Late-type  stars  (M  dwarfs)  and  brown  dwarfs 
also  exhibit  flaring  activity,  particularly  at  centimeter- 
and  meter  wavelengths  (;>,  300  MHz,  e.g.,  Ch'idell 
2002:  Hallinan  et  alJ[2007h.  Little  is  known  about 
the  emission  from  brown  dwarfs  at  LWDA  fre¬ 
quencies.  However,  the  most  sensitive  searches 
for  decameter-wavelength  emission  from  late-type 
stars  have,  at  best,  marginal  detections,  even  when 
those  stars  were  observed  to  be  flaring  at  shorter 
wavelengths  Ijackson  et  al.l  ( 1990h.  Even  these 
marginally  detected  flares  have  flux  densities  of 
order  1  Jy.  While  more  numerous  and  likely  to  be 
closer  than  solar- type  stars,  the  LWDA  sensitivity 
also  means  that  we  can  place  only  modest  limits 
on  the  likelihood  of  strong  flares  from  nearby  M 
dwarfs. 

Finally,  within  the  solar  system,  in  addition  to 
solar  radio  bursts,  Jupiter  also  produces  intense 
radio  bursts  due  to  an  electron  cyclotron  maser 
within  its  magnetosphere  ( Melrosel  20051.  This 
emission  cuts  off  above  around  40  MHz,  a  fre¬ 
quency  determined  by  the  ratio  between  the  cy¬ 
clotron  and  plasma  frequencies  as  a  function  of 
altitude  within  the  Jovian  magnetosphere. 


The  magnetospheric  emissions  from  solar  sys¬ 
tem  planets  and  the  discovery  of  extrasolar  plan¬ 
ets  have  motivated  a  number  of  both  theoretical 


(Zarka  et  al.  1997: 

Farrell  et  al.  1999t  Zarka  et  al. 

2001 

:  Farrell  et  al.1 

2004;  Lazio  et  al.  2004;  Stevens 

2005 

;  Griessmeier  et  al.  2005t  Zarka  20061  200tI 

Griessmeier  et  al. 

2007b 

ah  and  observational 

work  (Yantis  et  al.  1977 

Winglee  et  al.  1986| 

Bastian  et  al.  2000; lLazio  et  al.  20041:  Rvabov  et  al. 

2004 

:  George  &  Stevens 

2007t  Lazio  &  Farrell 

2007 


Smith  et  alJ2009ji~  on  magnetospheric  emis¬ 


sions  from  extrasolar  planets,  including  some  be¬ 
fore  the  confirmed  discovery  of  any  extrasolar 
planets.  Unfortunately,  even  the  most  optimistic 
predictions  for  extrasolar  planetary  radio  emission 
do  not  predict  flux  densities  in  excess  of  1  Jy  at  fre¬ 
quencies  near  100  MHz.  Moreover,  the  most  sen¬ 
sitive  searches  for  extrasolar  planetary  radio  emis¬ 
sion  near  74  MHz  place  sub- Jansky  constraints  on 
this  emission  (Lazio  fe  Farrelll  l2007t  lLazio  et  al 
200911. 


5.3.  Pulsar  Giant  Pulses 

We  consider  giant  pulses  from  radio  pulsars 
as  an  exemplar  of  a  short-duration,  intense  ra¬ 
dio  transient  that,  in  principle,  could  be  de- 
tected  throughout  the  Local  Group  by  the  LWDA 
(Gordes  et  al  l  2004a{l.  For  the  specific  case  of  the 


Rankin  et  al.  1970fl 


(Cornelia  et  al. 

1969 

Sallmen  et  al.  ( 

1999 

showed  that  the  pulses  are  indeed  broadband,  at 
least  over  the  frequency  range  610-4900  MHz. 

The  spectral  index  for  Crab  pulsar  giant  pulses 
is  steep  ( a  ~  —3.4,  Sv  oc  va),  and  current  ob¬ 
servations  find  no  upper  limit  to  the  amplitude  of 
giant  pulses  (Sallmen  et  al.  19991:  iLundgren  et  al.l 
1995ll .  However,  by  comparison  to  lower  frequency 
observations,  Lundgren  et  ah  (1199511  also  find  that 
the  rate  of  giant  pulses  decreases  with  decreas¬ 
ing  observational  frequency,  such  that  giant  pulses 
at  800  MHz  are  400  times  more  frequent  than 
at  146  MHz.  Scaling  both  the  rate  of  giant  pulses 
and  the  pulse  amplitudes  to  the  LWDA  frequency 
of  73.8  MHz,  we  estimate  that,  on  average,  the 
Crab  pulsar  will  produce  a  giant  pulse  once  every 
50  hr  with  a  flux  density  exceeding  1  MJy  (i.e., 
~  2000er  for  the  LWDA,  Table  HJ. 

At  higher  frequencies,  Crab  pulsar  giant  pulses 
have  extremely  narrow  widths,  often  being  unre- 
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solved  in  time.  With  a  pulse  period  of  approxi¬ 
mately  33  ms,  even  if  the  Crab  pulsar  or  a  pul¬ 
sar  like  it  produced  a  sufficiently  strong  pulse  to 
be  detected,  only  a  single  baseline  would  appear 
to  be  illuminated  by  our  data  acquisition  proce¬ 
dure  (51  ms  per  baseline,  !j3]).  However,  at  LWDA 
frequencies,  pulse  broadening  becomes  significant. 
Pulse  broadening  increases  the  duration  of  the 
pulse,  at  the  cost  of  decreasing  its  amplitude. 

In  addition  to  being  frequency  dependent, 
pulse  broadening  also  depends  upon  direction 
through  the  Galaxy.  Using  the  NE2001  model 
(Cordes  fc  Laziof2002ll.  we  estimate  that  the  mag¬ 
nitude  of  pulse  broadening  can  exceed  500  ms 
(~  10  LWDA  baselines  being  illuminated)  for  dis¬ 
tances  of  order  5  kpc  through  the  Galactic  disk, 
implying  Galactic  latitudes  |&|  <  5°.  Unfortu¬ 
nately,  these  low  Galactic  latitudes  constitute  only 
a  small  fraction  of  the  total  sky  (<  1%),  and  we 
also  blanked  the  Galactic  plane  as  part  of  our  ef¬ 
forts  to  remove  strong  sources  from  the  images  be¬ 
fore  conducting  the  statistical  tests.  Thus,  we  con¬ 
clude  that  while  giant-pulse  emitting  pulsars  seem 
capable  of  producing  sufficiently  strong  pulses  to 
be  detectable,  either  our  serial  data  acquisition 
procedure  would  have  not  been  sensitive  to  such 
pulses  or  significant  pulse  broadening  would  have 
rendered  them  too  faint  to  be  detectable. 


5.4.  Gamma-ray  Bursts  (GRBs) 


We  conclude  this  section  by  considering  radio 
pulses  from  GRBs,  as  an  exemplar  of  possible  ra- 
dio  pulses  of  extragalactic  origin.  IJsov  fc  Katz) 
(2000)  and  Sagiv  fc  Waxmanl  ( 20021  both  have 
predicted  that  GRBs  should  have  associated 
prompt  emission,  most  likely  below  100  MHz.  Pre¬ 
dicted  flux  densities  are  highly  model  dependent 
and  range  from  essentially  undetectable  to  in  ex¬ 
cess  of  1  MJy. 

Observationally,  there  have  been  a  number  of 
searches  for  radio  pulses  associated  with  GRBs. 
Cortiglioniet_al.  (1981),  Inzani  et  al.1  ( 19821).  and 
Amv  et  all  (119891  all  detected  some  dispersed  ra¬ 
dio  pulses,  but  found  no  convincing  associations 
with  GRBs.  iBalsanol  (19991)  found  a  dispersed  ra¬ 
dio  pulse  apparently  coincident  with  GRB  980329, 
but  it  was  narrowband,  which  has  led  to  it  be¬ 
ing  interpreted  as  due  to  terrestrial  interference. 
Various  searches  for  radio  pulses  associated  with 
GRBs  (including  precursor  pulses)  have  been  con¬ 


ducted  at  151  MHz 


Koranvi  et  al.  ( 1994.  19951): 


Dessenne  et  al.l  (119961).  Typical  upper  limits  have 


been  approximately  100  Jy. 


5.5.  Exotica 


One  of  the  motivations  for  this,  and  similar, 
searches  is  simply  that  previously  unknown  classes 
of  sources  may  be  discovered.  While  it  is  difficult 
to  assess  the  probability  of  an  unknown  class  of 
sources,  we  can  consider  extraterrestrial  transmit¬ 
ters  (Cocconi  fe  Morrisonlll959l)  as  an  example  of 
an  exotic  population  of  sources  o 

As  Loeb  fe  Zaldarriagal  (20071)  discuss,  much 
of  the  human-generated  radio  radiation  from  the 
Earth  is  emitted  in  the  range  50-400  MHz.  In¬ 
deed,  the  upper  end  of  the  LWDA’s  operational 
frequency  range  was  chosen  to  avoid  FM  ra¬ 
dio  broadcasts,  and  H4.ll  illustrates  the  detection 
of  TV  signals  with  the  LWDA.  In  contrast  to 


Loeb  fc  Zaldarriagal  (120071).  however,  we  consider 


short  duration  pulses,  such  as  might  originate  from 
a  transmitter  on  a  planet  rotating  into  or  out  of 
view,  whereas  they  consider  detecting  signals  in 
long  integrations  (~  1  month). 

The  LWDA’s  frequency  sub-bands  of  20  kHz  are 
well  matched  to  the  typical  bandwidths,  at  least 
for  terrestrial  transmitters  in  this  frequency  range. 
There  is  an  increasing  usage  of  so-called  spread 
spectrum  transmitters,  with  bandwidths  (much) 
larger  than  20  kHz.  However,  these  are  often  at 
higher  frequencies,  and  an  LWDA-like  instrument 
would  still  receive  all  of  the  power  within  its  re¬ 
ceived  bandwidth. 


Unfortunately,  with  a  500  Jy  sensitivity,  it  is  a 
relatively  simple  matter  to  show  that  a  transmit¬ 
ter  at  a  distance  of  10  pc  would  have  to  have  an 
effective  isotropic  radiated  power  (EIRP)  of  about 
1017  W  (170  dBW)  to  have  been  detectable  with 
the  LWDA.  By  contrast,  the  EIRP  of  even  pow¬ 
erful  military  radars,  such  as  the  Air  Force  Space 
Surveillance  System  (AFSSS),  only  approach  5  x 
1013  W  (137  dBW).  Further,  given  that  all  (or 
most)  of  the  stars  within  10  pc  are  cataloged,  it 
is  not  clear  that  an  all-sky  search,  such  as  the  one 
pursued  here,  is  the  most  effective  search  strat¬ 
egy,  unless  one  postulates  a  set  of  freely  floating 


4  Doing  so  is  also  appropriate  given  that  this  year  is  the 
50th  anniversary  of  the  first  search  for  radio  transmissions, 
F.  Drake’s  Project  OZMA. 
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“beacons.” 


5.6.  General  Considerations 


Examination  of  Figure  [8]  shows  that  an  LWDA- 
like  instrument  might  be  able  to  probe  a  significant 
region  in  this  phase  space,  and  previous  sections 
have  discussed  a  range  of  source  classes  that  the 
LWDA  might  have  been  able  to  detect,  in  princi¬ 
ple,  even  though  the  implementation  of  the  data 
acquisition  system  resulted  in  no  detections.  In 
addition  to  the  aforementioned  issue  that  a  more 
rapid  time  sampling  would  likely  be  required,  are 
there  other  considerations  for  a  dipole  phased  ar¬ 
ray? 


Much  of  the  phase  space  accessible,  or  poten¬ 
tially  accessible,  to  an  LWDA-like  instrument 
is  for  coherent  transients,  nominally  taken  to 
be  those  with  brightness  temperatures  in  excess 
of  1012  K.  Such  transients  will  necessarily  evolve 
rapidly.  ICordes  fe  McLauehlinl  ( 2003h  have  dis¬ 
cussed  searches  for  fast  radio  transients.  A  key 
consideration  for  such  searches  is  the  effective  time 
resolution,  which  is  determined  not  only  by  the  in¬ 
strumental  characteristics  (receiver  bandpass)  but 
also  by  dispersion  smearing  and  radio-wave  scat¬ 
tering. 

Following  ICordes  &  McLaughlin  (l2003lh  we 
consider  two  time  scales  related  to  the  time  reso¬ 
lution,  and  implied  processing,  for  searches  with 
an  LWDA-like  instrument.  The  highest  time  reso¬ 
lution  that  the  instrument  would  be  able  to  obtain 
is  simply 

A tAv  ~  (Ai/)"1,  (3) 


for  a  receiver  frequency  channelization  or  sub¬ 
band  of  Ai/,  which  for  the  LWDA  (Au  =  20  kHz) 
was  50  /is.  Because  of  the  dispersion  relation  of 
the  interstellar  medium,  probing  to  any  distance 
for  fast  pulses  requires  de-dispersion.  As  an  illus¬ 
tration,  we  consider  how  far  one  could  probe,  ig¬ 
noring  dispersion.  We  take  as  a  criterion  that  the 
dispersion  across  the  receiver  sub-band  produces  a 
smearing  comparable  to  time  resolution  implied  by 
Equation  0.  Then,  from  Cordes  &  McLaughlinl 
(l2003h. 


(5DM  =  120  pc  cm  3  (Ai^Hz)  24hz.  (4) 

where  A^kHz  is  the  receiver  sub-band  in  units  of 
kHz,  and  ^mhz  is  the  central  frequency  in  units 


of  MHz,  and  the  error  in  the  dispersion  mea¬ 
sure  (5DM  is  in  the  canonical  units  of  pc  cm'3. 
Alternately,  <5DM  can  be  considered  the  incre¬ 
ment  in  dispersion  measure  (DM)  for  searches 
that  incorporate  searching  through  dispersion. 
For  the  LWDA,  (5DM  =  0.1  pc  cm-3,  imply¬ 
ing  that  a  search  for  fast  transients  with  the 
LWDA  that  did  not  incorporate  de-dispersion 
could  only  probe  to  an  effective  DM  of  this  value. 
For  reference,  PSR  J0108— 1431  has  the  lowest 
known  DM,  2.38  pc  cm~3,  at  a  distance  of  240  pc 
(Manchester  et  al.  2005:  Deller  et  al.1  2009li. 


Conducting  a  search  that  incorporated  compen¬ 
sation  for  dispersion  smearing  would  effectively 
mean  conducting  a  de-dispersion  search  on  each 
resolution  element  (beam)  within  the  held  of  view 
(Cordesl  12009).  In  contrast  to  proposed  searches 
at  higher  frequencies  (~  1  GHz),  the  much  lower 
resolution  at  LWDA  frequencies  makes  this  a  less 
onerous  task.  For  instance,  with  the  LWDA’s  res¬ 
olution  of  12°,  there  are  only  about  150  indepen¬ 
dent  resolution  elements  in  the  entire  sky.  Further, 
searching  each  is  of  course  an  embarrassingly  par¬ 
allel  problem. 


6.  Conclusions 

We  have  described  the  Long  Wavelength  Demon¬ 
strator  Array  (LWDA)  and  its  operation  as  an 
all-sky  transient  monitor.  A  16-element  dipole 
phased  array,  operating  over  the  frequency  range 
60-80  MHz,  the  LWDA  was  used  as  a  technical 
test  bed  for  the  LWA.  We  have  described  the  sig¬ 
nal  how  from  the  antennas  to  the  receivers  and  to 
a  simple  software  correlator  as  well  as  the  config¬ 
uration  of  the  dipoles. 

The  individual  dipoles  of  the  LWDA  had  a  held 
of  view  comparable  to  the  entire  sky.  In  late  2006 
and  early  2007,  we  used  the  LWDA  as  a  tran¬ 
sient  search  instrument  by  making  a  series  of  all¬ 
sky  images  at  frequencies  of  61  and  73.8  MHz. 
The  61  MHz  observations  were  designed  explicitly 
to  search  for  reflections  from  ionized  meteor  trail 
during  the  Leonid  meteor  shower  in  2006.  The 
73.8  MHz  observations  were  designed  to  search  for 
astronomical  transients,  and  we  acquired  a  total 
of  106  hr  of  data,  with  a  time  sampling  ranging 
from  2  to  5  minutes  between  images.  We  were  able 
to  detect  solar  hares  and,  by  utilizing  a  special- 
purpose  mode,  TV  reflections  off  ionized  meteor 
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trails  during  the  2006  Leonid  meteor  shower.  We 
detected  no  radio  transients  outside  of  the  so¬ 
lar  system  above  a  flux  density  limit  of  500  Jy, 
equivalent  to  a  limit  of  no  more  than  about  10-2 
events  yr-1  deg-2,  having  a  pulse  energy  density 
>  1.5  x  10-20  J  m-2  Hz-1  at  73.8  MHz  for  pulse 
widths  of  about  300  s.  This  event  rate  is  compa¬ 
rable  to  many  existing  limits  from  previous  all-sky 
surveys,  but  at  a  lower  frequency  than  most  pre¬ 
vious  all-sky  searches. 

There  are  a  number  of  emerging  arrays  (e.g., 
MWA,  LOFAR,  LWA)  for  which  all-sky,  or  at  least 
very  wide-held,  imaging  will  be  a  capability.  All  of 
these  arrays  will  offer  at  least  two  advantages  with 
respect  to  the  LWDA.  First,  all  are  anticipated 
to  be  operational  for  a  far  longer  time  than  was 
the  LWDA,  which  should  improve  upon  the  event 
rates  substantially.  Second,  all  of  these  arrays  will 
have  a  higher  angular  resolution  imaging  capabil¬ 
ity,  so  that  if  a  transient  were  to  be  detected,  much 
higher  precision  astrometric  information  could  be 
obtained. 

The  LWDA  was  a  joint  project  of  NRL,  the 
Applied  Research  Laboratories  of  the  University 
of  Texas  at  Austin,  and  the  University  of  New 
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with  the  infrastructure,  construction,  and  main¬ 
tenance  of  the  LWDA  and  the  many  staff  mem¬ 
bers  of  the  NRAO  who  provided  technical  sup¬ 
port  in  establishing  the  LWDA  near  the  VLA.  We 
thank  K.  Weiler,  D.  Munton,  and  L.  J  Rickard 
for  their  guidance  during  various  stages  of  the 
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The  LUNAR  consortium,  headquartered  at  the 
University  of  Colorado,  is  funded  by  the  NASA 
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